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Abstract. The celebration of AUan Kaufman's 80th birthday was an occasion to reflect on 
a career that has stimulated the mutual exchange of ideas (or memes in the terminology of 
Richard Dawkins) between many researchers. This paper will revisit a meme AUan encountered 
in his early career in magnetohydrodynamics, the continuation of a magnetohydrodynamic mode 
through a singularity, and will also mention other problems where AUan's work has had a 
powerful cross-fertilizing effect in plasma physics and other areas of physics and mathematics. 



1. Introduction 

Richard Dawkins [B 12]) in discussing the evolution of ideas, introduced the concept of memes 
in the foUowing terms ". . . new rephcators are not DNA . . . patterns of information that can 
thrive only in brains . . . or books, computers and so on ... cahed memes to distinguish them 
from genes . . . passed from brain to brain ... As they propagate they can change mutate . . . 
memic evolution ..." 

Applying and extending this biological analogy to the evolution of science, we can regard 
scientific papers as the "organisms" in which multiple memes are expressed, and say that cross- 
fertilization leads to hybrid vigour. 

One of Allan Kaufman's greatest contributions to plasma physics has been his role in adapting 
and developing powerful theoretical and mathematical methods, applying them to plasma 
physics problems, and propagating these ideas to foUowing generations through his ability as a 
great teacher. 

My (RLD) research interests have intersected, and continue to intersect, with Allan's varied 
research interests. In this paper I have chosen to illustrate this with a few examples from 
magnetohydrodynamics (MHD). This may seem a little unexpected, because Allan is not usually 
thought of as an MHD theorist, but in section[2]l revisit a problem he was familiar with in his 
early career, namely the continuation (or lack thereof ) of solutions of the Newcomb equation for 
linearized ideal marginal MHD modes. The disconnection between solutions on either side of a 
mode rational surface is a tricky point to explain to newcomers in the field and I present a new 
approach which I hope makes this phenomenon clearer. 

In section[3]l mention some other MHD problems I have encountered which have intersections 
with Allan's research. 



2. The connection problem at a singular point of the Newcomb equation 

It is probably not widely known that some of Allan's earliest work was on magnetohydrodynamic 
(MHD) stability theory [3 . Newcomb's famous paper on ideal-MHD stability theory in 
cylindrical geometry [4 acknowledges Allan's critical reading of the manuscript while his Ref. 24 
is the note "M. N. ROSENBLUTH (private communication through A. N. Kaufman)," showing 
that Allan's cross-pollinating role in plasma physics started very early on. 

It happens that these papers are highly relevant to our current research |5l [6l [7| on equilibrium 
and stability in a multiple-region relaxed-MHD plasma model. This is motivated by the desire to 
construct a well-posed three-dimensional MHD equilibrium theory, based on Taylor relaxation in 
regions separated by arbitrarily thin ideal-MHD toroidal surfaces that act as barriers to field-line 
chaos. Not only does this involve the Hamiltonian nonlinear dynamics of the field-line flow, but 
the Hamiltonian of the field lines must itself be determined self-consistently using MHD theory. 

To understand the nature of these ideal-MHD barriers better we have recently |8l [9] returned 
to the cylindrical limit and looked at the problem of whether the assumed barrier between two 
relaxed regions of different pressure can be constructed as the zero-width limit of a finite-width 
ideal-MHD region with a suitably chosen physical pressure profile. [The criterion for physicality 
of the pressure profile is that the pressure P(r) must be non-negative, where r is the distance 
from the 2;-axis.] This problem was first studied by Newcomb and Kaufman [3j ;4J. 

The "Newcomb equation" |5j is satisfied at marginal stability (i. e. when the growth rate 
7 is zero) by ^(r) = ^-e^, where $, oc expz(m^ -|- kz) is the plasma displacement away from 
equilibrium, 9 being the angle about the 2:-axis and e,, the unit vector in the radial direction: 

where 
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For given integers m and n = —Rk 7^ O, a mode rational surface is one where q{rs) = m/n, 
with q{r) = rBz/RBg and R being the periodicity length of the cylinder. At such a rational 
surface r^, ([T| has a singular point because the coefficient of the highest derivative, f (r), 
vanishes there. Defining x = r — r^, the Taylor expansions of / and g about r g are of the 
form /(r) = f"{rs)x^/2 + 0{x^) and g{r) = g{rs) + 0(x), with f "{r s) > O provided q'{rs) / O 
and g {r s) / O provided P' {r g) 7^ 0. The general solution of ([T| on intervals not including r s is 
found as a linear superposition of two functions whose Frobenius expansions (see e. g. |10| ITT] ) 
begin with the fractional powers where 
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We assume /i is real, otherwise the system is interchange-unstable by the Suydam criterion. 
The Solutions with leading terms and were called by Newcomb the large and 

small Solutions, respectively, the large solution being dominant as x — > 0. 

In this paper we discuss a point raised by Newcomb jH p. 241]: "In general it is not possible 
to continue an Euler-Lagrange solution past a singular point." To this he adds the footnote: 



"This is true only because is a real variable. If it were complex we could go around the singular 
point by analytic continuation, but the resulting solutions would generally be multivalued." The 
issue is interesting in the context of the present paper because this meme is perhaps expressed 
in a mutated form in Allan's later papers on mode conversion (e. g. [E]), which also involves 
continuation through singularities. 

The natural framework in which to understand the singular solutions of the Newcomb 
equation is generalized function theory because the solutions are not defined pointwise (being 
undefined at r = r^). In our previous work [13l [Hl [15], using the singular solutions of the 
Newcomb equation or the two-dimensional generalization of it due to Bineau [16] , we have 
followed the approach and notation of of Gel'fand and Shilov [17] based on inner products with 
all members of a space of smooth test functions. 

However this approach is rather abstract and difficult to visualize, so in the current paper we 
present an alternative approach, where the required generalized functions are defined as limits of 
sequences, parametrized by 5, of smooth functions that tend to weak solutions of the Newcomb 
equation as 5 — > +0. This is in the spirit of Lighthill's [18] approach to generalized function 
theory. 

To do this we study a model Newcomb equation, in which /(x) = x^ and g{x) = — 1/4. 
The weak solution space is spanned by the four Gel'fand-Shilov generalized functions 
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where sgnx denotes the sign of x. On the left-hand sides we have used the notation of Dewar 
and Persson pT| in which the large and small Frobenius solutions are indicated by using upper 
and lower case greek, respectively; odd and even parity solutions are indicated by subscript — 
and +, respectively. 

Alternatively, we may use the "single-sided" generalized functions = ^d^l"^ it |x|'^sgnx) 
whose support is the positive (x > 0) or negative (x < 0) half line. When restricted to their 
domains of support these single-sided solutions correspond to the classical pointwise solutions 
Newcomb meant when he said solutions of the Newcomb equation could not be continued through 
the singular point. As the domains of support are to the left (L) and right (R) of the origin, we 
denote the corresponding generalized function solutions of the Newcomb by subscripts L and R 

E^{x) ^ 

Cl^{x) ^ \x\-'/'^^. (6) 

In ideal MHD the large solutions are rejected as they are not square integrable (hence give 
infinite kinetic energy) , thus reducing the dimensionality of the solution space to 2 as expected 
for second-order differential equations. (However, the large solutions are essential in resistive 
stability theory for matching asymptotically to the resistive internal layer in the neighbourhood 
of the singular point.) 

To apply the Lighthill approach we use a regularization appropriate to ideal-MHD stability 
studies, in which the only regularizing effect is that of inertia from the mass density p. This 
comes into play when the growth rate 7 is nonzero, adding to / a positive term proportional to 



/97^ (see e. g. [19 ) and thus removing the singularity but approximating the Newcomb equation 
arbitrarily closely as 7 ^ 0. Thus we study the regularized model Newcomb eguation 



whose general solution is a hnear combination of the respectively odd and even functions 
Im P^_i/2{ix/6) and ReP^_i/2{ix/6), where Pu{z) = Pl^{z) denotes a Legendre function of 
the first kind [20, chapter 8]. 
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Figure 1. Left panel: regularized odd solution ^^''(x|0.001) (solid line) and unregularized odd 
solution C^{x) (dashed line). Right panel: regularized even solution ^^''(x|0.005) (solid line) 
and unregularized even solution £,^{x) (dashed line). In both cases /u = 0.4. 



By expanding Pi^-i/2{ix/S) in powers of 6 and choosing the superposition coefficients 
appropriately we can now define Lighthill approximating functions for the odd and even small 
Solutions 

' ^ - sin(f (;.-i))r(2^)'"^'^"-H5;' 



+^ ' ^ cos^ /i-i r(2/.) ^-^v^; 
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Figure [T] shows examples of the unregularized solutions defined in ^ and the regularized 
Solutions defined in ([s]) and ([9|. 

The leading terms of the asymptotic expansions of these two functions in powers of are 
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where the factors c±(^) are defined by 
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c (a) ^ o2. r(/i+i)^r(-2/x)sin(f (/.+ !)) 
r(i-/.)V(2^) sin(f (/.-i)) 
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The 0(1) terms are the odd and even generalized function small solutions as required. The 
O {6'^'^) terms involve the large solutions of the same parity, with coefficients that vanish as 
5 +0 because our regularization is based on ideal MHD. 

The special case /i ^ 1/2 is important because it corresponds to the zero-/? limit, or flattening 
pressure at a rational surface at arbitrary /3. In this limit we find 

ax\S) = sgnx - —x-^ + 0{6^), 
vr 

^f{x\S) = l + 0(<52). (12) 
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Figure 2. Left panel: regularized left solution ^^'^(slO.OOl) (solid line) and unregularized left 
solution Ci!^{x) (dashed line). Right panel: regularized right solution ^^''(2;|0.001) (solid line) 
and unregularized right solution ^^"^(ic) (dashed line). In both cases /j, = 0.4. 



We are now ready to regularize Newcomb's "disconnected" solutions by combining the 
regularized odd and even parity solutions in an analogous way to that used for unregularized 
solutions in (|6|, 
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so that the 0{6^) terms cancel to the right or left of the origin for ^i]'^{x\S) and ^^'^(j;|(5), 
respectively. 

Figure |2] shows examples of the unregularized solutions defined in ([6| and the regularized 
solutions defined in (13) and (14). It is seen that, for any finite 5, the solutions do connect 

ITol) 



across the singular point at the origin. By using ( |10[ ) in (13) and (14) we see that the L and 
R solutions decay rapidly (like the large solution) as x ^ ±00 with an amplitude that tends to 
zero as (5 — > +0, so that their support becomes the negative or positive half line in the limit. 



3. Other MHD intersections 

The most obvious point of intersection between AUan Kaufman's and my (RLD) research 
interests is in the oscillation centre theory and the theory of wave action, which goes back 
to an MHD paper |21j . This has developed in various ways that other contributors to this 
Proceedings will no doubt touch on in much more detail. However, sticking to the MHD theme, 
it could be said that the theory of flux-minimizing surfaces |22| l23l [2^ [25l [26] is a mutant form of 



oscillation centre theory, and this has close links with our current research on three-dimensional 
MHD equilibrium theory [H] . 

Another intersecton concerns the theory of quantum chaos, a field in which Ahan made 
an important contribution [27 , which is better known outside the field of plasma physics than 
within it. Indeed there has been very little in the plasma physics literature about quantum chaos, 
perhaps because of the name. However, this field is really about WKB theory for arbitrary waves 
in the case of non-integrable ray dynamics, and in a number of papers over the past few years 
we have applied methods from this field to the study of the ideal-MHD spectrum in a stellarator 



4. Conclusion 

In this brief paper we have revisited a basic problem in the the theory of ideal MHD stability 
theory that has connections with Allan Kaufman's early research and our current research. In 
his long career Allan has been a role model and inspiration to many, and particularly to the first 
author of this paper. 
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